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Abstract: Following up on recent work in the context of ordinary fluids, we study 
the equilibrium partition function of a 3+1 dimensional superfluid on an arbitrary 
stationary background spacetime, and with arbitrary stationary background gauge 
fields, in the long wavelength expansion. We argue that this partition function is 
generated by a 3 dimensional Euclidean effective action for the massless Goldstone 
field. We parameterize the general form of this action at first order in the derivative 
expansion. We demonstrate that the constitutive relations of relativistic superfluid 
hydrodynamics are significantly constrained by the requirement of consistency with 
such an effective action. At first order in the derivative expansion we demonstrate that 
the resultant constraints on constitutive relations coincide precisely with the equalities 
between hydrodynamical transport coefficients recently derived from the second law of 
thermodynamics. 
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1. Introduction 

Hydrodynamics is the long distance effective description of locally thermalized systems. 
The variables of hydrodynamics are local values of the temperature, chemical poten- 
tial, velocity and other relevant thermodynamical order parameters. The equations 
of hydrodynamics are the universal laws of conservation of the stress tensor and the 
charge current. These equations may be used to describe the propagation of a fluid 
in an arbitrary weakly curved background metric, and with arbitrary slowly varying 



background gauge fields. Within the hydrodynamical approximation, the stress tensor 
and charge current of the fluid are expressed as functions of thermodynamical and 
background fields; the formulas through which this is achieved are referred to as the 
constitutive relations of the hydrodynamical system. 

The hydrodynamical constitutive relations of any given system may, in principle, 
be determined by a detailed study of the dynamics of the theory. For strongly coupled 
quantum field theories, however, the required calculations usually cannot be practically 
executed. 1 Given this state of affairs, it is clearly of interest to have a complete 
parameterization of the most general hydrodynamical constitutive relations allowed on 
general grounds. Such a characterization would constitute a satisfactory framework 
for the the theory classical hydrodynamics viewed as an autonomous long wavelength 
effective theory. 

The constitutive relations of relativistic hydrodynamics are specified in an expan- 
sion in derivatives of the local thermodynamical fields and background fields. At any 
order in this expansion, Lorentz invariance determines the constitutive relations up to 
a finite number of functions of the scalar thermodynamical fields (e.g. local tempera- 
ture and chemical potential). It turns out, however, that other considerations further 
constrain the constitutive relations. These constraints are of two sorts. The first, and 
more important, set of constraints asserts relations between the apparently indepen- 
dent functions that appear in the most general Lorentz allowed constitutive relations. 
Such constraints cut down the number of free functions in the equations of hydrody- 
namics (at any order in the derivative expansion). We refer to constraints of this first 
sort as 'equality type' constraints; they are the primary focus of the current paper. A 
second, milder form of constraints assert inequalities for the free functions that appear 
in constitutive relations. These constraints do not reduce the number of free functions 
in constitutive relations but merely bound these functions. In this paper we will have 
nothing to say about this second class of constraints. 

One method for obtaining constraints on constitutive relations was outlined in 
the classic text book of Landau and Lifshitz || and refined in later studies ( see e.g. 
H H) 0> H> H3 f° r recen t work). This method is based on the assumption that 
consistent equations of hydrodynamics come equipped with an entropy current. Like 
the conserved currents, the entropy current is a function of the local thermodynamical 
fields. The key dynamical assumption is that the divergence of this entropy current 
is point wise (in spacetime) positive semi definite for every conceivable fluid flow. We 
refer to the existence of such a positive divergence entropy current - without making 

1 The only exceptions that we are aware of lie within the fluid gravity map ( jjj], see see ^ ||] for 
reviews) of the AdS/CFT correspondence of string theory. The constitutive relations of field theories 
with a known dual description are rather easily determined following the procedure first described in 



any prior assumptions as to its functional form - as a local version of the second law of 
thermodynamics. The local second law clearly guarantees that net thermodynamical 
entropy increases in any fluid flow that starts and ends in equilibrium. The requirement 
that entropy increase for fluid flows perturbed by arbitrary local sources suggests that 
the local form of the second law is also a necessary consequence of the second law. 2 

At low orders in the derivative expansion it has been demonstrated that the lo- 
cal version of the second law of thermodynamics yields powerful constraints (of both 
the equality and inequality sort) on otherwise unrelated free functions in constitutive 
relations. The detailed form of these constraints has been worked out for uncharged 
relativistic fluids up to second order in the derivative expansion |10] and for (parity 
non preserving) charged fluids to first order in the derivative expansion (see |5], [11], |6j 
for 3+1 dimensional fluids and || for 2+1 dimensional fluids). 

Very recently, a second systematic method for constraining the constitutive re- 
lations for fluids was described in |TJ1 [13]] 3 . These constraints follow from the very 
reasonable demand that the hydrodynamical equations must always admit equilibrium 
solutions for arbitrary stationary background field configurations, and moreover that 
the values of conserved charges in equilibrium must be consistent with the existence of 
an equilibrium partition function. For both uncharged fluids at second order as well 
as charged parity non invariant (and potentially anomalous) fluids at first order, the 
equality type constraints obtained from this method agree, in full detail [|12|, [13|] with 
those obtained from the local form of the second law of thermodynamics described 
above. It has been conjectured []12 that this agreement persists to all orders in the 
derivative expansion; however there is as yet no proof of this conjecture. 

The equations of charged hydrodynamics are modified when the charge symmetry 
of the system is spontaneously broken by the condensation of a charged operator in 
thermal equilibrium. The effective description of such systems has new hydrodynamical 
degrees of freedom whose origin lies in the Goldstone mode of the charge condensate. 
The resultant hydrodynamical equations are referred to as the equations of superfluid 
hydrodynamics, and are the subject of the current paper. 

More particularly in this paper we study 's' wave superfluid hydrodynamics, i.e. 
the hydrodynamics of a system whose charge condensate is a complex scalar operator. 
We study the constraints on the equations of first order 's' wave superfluid hydro- 
dynamics imposed by the requirement that these equations admit equilibrium under 



2 It is at least suggestive that the fluid dynamics generated by the fluid gravity map of the AdS / CFT 
correspondence does indeed always obey the local form of the second law, at least when the dual bulk 
description is given by the equations of two derivative Einstein gravity with matter that obeys the 
null energy conditions. 

3 See also |l4|, for earlier related work, [^6| for closely related discussions and |l7], Q for follow 
up work. 



appropriate situations, and that the charge currents in equilibrium agree with those 
from an appropriate partition function. We do not assume that the superfluids we 
study necessarily preserve either parity or time reversal invariance. 

As we explain in section below, the general analysis presented in this paper 
closely follows that of fl2| (for the case of ordinary, i.e. not 'super' fluids) with one 
important difference. The Euclidean partition function for a superfluid in an arbitrary 
background 4 is determined by an effective field theory that includes a massless mode: 
the Goldstone boson of the theory. This effective field theory is local, and may usefully 
be studied in the derivative expansion. However the partition function that follows after 
integrating out the Goldstone boson is neither local nor simple. As we explain below, 
the study of the local effective action of the Goldstone boson (rather than the partition 
function itself) allows us to usefully constrain the constitutive relations of superfluid 
hydrodynamics. In this paper we present a careful derivation of the relations between 
otherwise independent transport functions that follow from such a study. 

Constraints on the constitutive relations of first order superfluid hydrodynamics 
have previously been obtained using the local form of the second law in |2(| |7], |[ || 
for the case of time reversal invariant superfluids. In this paper we generalize the 
derivation of || to include the study of superfluids that do not preserve time reversal 
invariance. We then compare the results obtained from the two different methods; i.e. 
the constraints that follow from the requirement of existence of equilibrium and those 
that follow from the local second law. As in the case of ordinary (i.e. non super) 
fluids we find perfect agreement between the equality type constraints obtained from 
these two apparently distinct methods. Our results supply further evidence for the 
conjecture that the equality type constraints from these two methods agree in a wide 
range of hydrodynamical contexts and to all orders in the derivative expansion. A proof 
of this conjecture would go some way towards proving the local form of the second law, 
and would permit the demystification of this law in a hydrodynamical context. 

While the work reported in this paper is purely hydrodynamical and nowhere uses 
AdS/CFT, much of the motivation for this work lies within the fluid gravity map of 
AdS/CFT. The status of the second law of thermodynamics for theories of gravity that 
include higher derivative corrections to the Einstein Lagrangian remains unclear. In 
particular it has never been proved that the Hawking area increase theorem general- 
izes to a Wald entropy increase theorem for arbitrary higher derivative corrections to 
Einstein's gravity. If the interplay between the existence of equilibrium in appropriate 
circumstances and entropy increase can be proved on general grounds in a hydrody- 
namical context, then it seems likely that the lessons learnt can be taken over to the 
study of entropy increase in higher derivative gravity (at least for asymptotically AdS 



4 See [|l3[ ^9) for a discussion of this partition function at the perfect fluid level. 



space) via the fluid gravity map. This could lead to a proof of a Wald entropy in- 
crease theorem under appropriate conditions on the higher derivative corrections of 
the gravitational system. 



2. Equilibrium effective action for the Goldstone mode 
2.1 The question addressed 

In this section we study an s wave superfluid propagating on the stationary background 
metric 

ds 2 = G^dx^dx" = -e 2a{s) (dt + a; (f) ate*) 2 + g ij (x)dx i dx j (2.1) 
and background gauge field 

A = A (x)dx° + A i (x)dx i (2.2) 
Below we will often work in terms of the modified gauge fields 



Ai = Ai — A ai 
A = Ao + jj 



(2.3) 



All background fields above are assumed to vary slowly; we work in an expansion in 
derivatives of these fields. We address the following question: what is the most general 
allowed form of the partition function 

Z = Tre T o (2.4) 

as a function of the background fields a, ai, gij, Ao and Ai in a systematic derivative 
expansion? 

2.2 The partition function for charged (non super) fluids 

The analogous question was studied for the case of an ordinary (non super) charged 



fluid in |12l . It was demonstrated that to first order in the derivative expansion the 
most general allowed form of the partition function for an ordinary charged fluid on 
the background ( |2~T|) , (|2.2|) is given by 

W = \nZ = W° + WL + W 1 anom 

V9^(Toe- ff , e -M„) 
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where P(T, (j) is the thermodynamical pressure of the system as a function of its 
temperature and chemical potential and Co, C\ C2 and C are all constants. The 
constant C specifies the covariant U(l) 3 anomaly via the equation 

9/ = 4*(fAF) (2.6) 

o 

The constants Co, C\ and C2 do not (yet) have similar interpretations. It was demon- 
strated that Co = Ci = in any system that respects CPT invariance. 

Notice that the result Q2.5|) for the partition function of an ordinary (non super) 
fluid is a local function of the background sources gij, aj, a, A and A4. Locality 
is a direct consequence of the fact that the path integral that computes the partition 
function ( |2.4|) has a unique hydrodynamical saddle point (as opposed to a moduli space 
of saddle points). As a consequence the partition function is generically 5 computed 
by a path integral over an action with no massless fields. It follows that the result is 
local on length scales large compared to the inverse mass gap in the action (this mass 
gap is sometimes referred to as a static screening length of the 4 d thermal system) 6 . 

2.3 Euclidean action for the Goldstone mode for superfluids 

Unlike an ordinary charged fluid, the equilibrium configuration of a superfluid in the 
background Q2.1|) is not unique. As superfluids break the global 27" (1) symmetry, every 
background admits at least a one parameter set of equilibrium configurations that differ 
by a constant shift in the phase of the expectation value of the condensed scalar. It 
follows that the path integral that computes ( |2.4|) has a zero mode (the phase of the 
scalar condensate). Consequently, the partition function fl2.4|) , is not a local function 
of the background source fields. Instead this partition function is generated by a local 
three dimensional field theory of the dynamical phase field (j). 

The dynamics of the Goldstone boson in general, governed by a 3d massless quan- 
tum field theory. In this paper, however, we focus on field theories in an appropriate 
large N limit (such as theories with matrix degrees of freedom in the t' Hooft limit). 
In such a limit the effective action for the Goldstone boson is multiplied by a suitable 
positive power of N (the factor is iV 2 in the t'Hooft limit mentioned above). As a 
consequence Goldstone dynamics is effectively classical in the large N limit. Quantum 
corrections to this classical answer, which are suppressed by appropriate powers of N 
(this power is in the t'Hooft limit), may have very interesting structure, see e.g. 
for related work. We leave their study to future work. 7 

5 Non hydrodynamical massless modes occur when the system is tuned to a second order phase 
transition. We assume in what follows that our system has not been tuned to such a phase transition. 
We leave the study of this interesting special case to future work. 

6 We thank K. Jensen for discussions on this topic 

7 We thank K. Jensen for discussions on this topic. 
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In principle, the partition function ( |2.4j ) for the superfluid may be obtained from 
the corresponding local effective action by integrating out the Goldstone boson (i.e. 
solving its equation of motion and plugging the solution back into the action). 8 In 
practice the implementation of this procedure requires the solution of a nonlinear par- 
tial differential equation. Moreover, even if one could solve this equation the resultant 
partition function would be highly nonlocal. A direct analysis of the partition func- 
tion itself seems neither easy nor particularly useful. In order to obtain constraints on 
the equations of superfluid hydrodynamics below we will work directly with the local 
effective action for the Goldstone mode rather than the final result for the partition 
function. 

The requirements of gauge invariance significantly constrain the form of Goldstone 
effective action. Let <fi denote the phase of the scalar condensate. Under a gauge 
transformation Ai — > A% + dia, 4> transforms as <fi + a - It follows that the effective 
action can only depend on the combination 

ii = -84 + Ai 

Note that £ M like A^, is a field of zero order in the derivative expansion 9 . 

The local field theory for the Goldstone boson must also enjoy invariance under 
Kaluza Klein gauge transformations (a^ — > Oj — cVy, see subsection 2.2 of JL2] for 
details). For this reason we work with the Kaluza Klein invariant fields 



We also define 



and define 



Ci = Zi-aiA Q = -84 + A. (2.7) 



X = T = -U" = & - 9 V (i(r (2-8) 
2.4 The Goldstone action for perfect superfluid hydrodynamics 

As we have explained above, the euclidean partition function for our system is generated 
by an effective action S for the Goldstone field 0. This Goldstone action may be 
expanded in a power series in derivatives. 

S = S + S 1 + S 2 . . . (2.9) 



8 If the Euclidean 3 dimensional manifold we work on is compact and we demand single valuedness 
of the field 4> then it is plausible that the solution to the <j> equation of motion is (at least generically) 
unique, see below. 

9 This means that the phase field <j> is of — 1 order in derivatives; this observation does not invali- 
date the derivative expansion as gauge invariant physical quantities are functions only of £ M and not 
independently of 0. 



At lowest (zero) order in the derivative expansion symmetries constrain the Goldstone 
boson effective action to take the form 10 



S = J d 3 Xy/g^P(f,fi,x)- 



T = T oe " (2.10) 
fi = A Q e^ 

vf 1 = (l,0,0,0)e" a 

where P is an arbitrary function whose thermodynamical significance we will soon 
discover, and x was defined in ( |2.8| ). The fields T, fi and are the values of the 
hydrodynamical temperature, chemical potential and velocity fields in equilibrium at 
zeroth order in the derivative expansion (see [12]). 

In the classical (or large N) limit adopted throughout this paper, the partition 
function Z of our system is obtained by evaluating the Goldstone action on shell. Let 
the solution to the equation of motion be denoted by 

G(s) = £*(*)• 

Then the partition function is given by 

\nZ = S(Cr(x)) (2.11) 

At lowest order in the derivative expansion, the action fl2.10|) depends only on first 
derivatives of the massless field <f>. Varying this action w.r.t. 

f e a BP 

J T dx 

= - j fx^-djiV^GrfW (2.12) 



yields 



where 



djiV^Gft 1 ) = vW(/e) = V, UA = 0. (2.13) 



dx 

Note this equation of motion is of second order in derivatives of the field <fi. 11 Plugging 
the solution to ( j2.13|) back into the (p.12 ) in principle yields an explicit though com- 



10 r 



The action ( |2. 10|) was already presented in J_3[ . The presentation of this subsection differs from 



|13| only in the emphasis that <f> be regarded as a dynamical held in (2.10), rather than a background 
like f. For related discussions on effective action for superfluid, see for example |26|, |l9|| . 

1:L The formal similarity of ( |2.13| ) to the equation V 2 = (where the Laplacian is taken in an 



appropriately rescaled metric) suggests that (2.12) has a unique solution on a compact manifold (up 



to constant shift in (j)) provided that (f> is required to be single valued and smooth on this manifold. 
However we do not have a proof of this statement. 



plicated and nonlocal expression for the partition function of the system as a function 
of source fields. 

The stress tensor and charge current that follow from the action fl2.10|) may be 



computed in a straightforward manner using the formulas listed in eqs.(2.16) of ||12|| ; 
they are given by 
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The gauge and diffeomorphism invariance of the action (|2.10| ) ensure the stress tensor 
and charge current described above are automatically conserved onshell (i.e. upon 
imposing the equation of motion ( |2.13| )). 

The complicated looking expressions for the conserved currents ( |2.14| ) may actually 
be summarized in a remarkably simple form as 



= (e + P)u^u u + Pg» v + f^C 



(2.15) 



where u was defined in ( p.lOj ) and all terms on the RHS of ( |2.15[ ) are evaluated on the 
zero order equilibrium solutions T(x) = T and n(x) = ft, defined in ( |2.10| ) and the 
functions e, s and q are defined in terms of the pressure p by the equations 



e + P = sT + qfi 

dP = sdT + qdfi + ^fdx 



(2.16) 



( |2.15| ) and (|3.2| ) are precisely the Landau-Tisza constitutive relations of superfluid 
hydrodynamics. 



2.5 The Goldstone Action at first order in derivatives 



One derivative corrections to the Goldstone action ( |2.10| ) may be divided into parity 
even and parity odd terms. We consider these in turn. 



2.5.1 Parity even one derivative corrections 

The most general parity preserving one derivative correction to ( [2. ICQ is given by 



n c I oevet 

ST en = I d 3 yVg 



where T was defined in (2.10), 



4(c.d)f +|(c.d)*>-/ 3 v i Uc 



A A 

v= t = % 



(2.17) 



and 



f^f^f^X 2 ) (2 = 1. -.3) 
are arbitrary functions while / was defined in the previous subsection 

dP 



f(T,u J C 2 ) = -2 



Two remarks are in order 



1. In ( [2.1 7|) the unspecified function f 3 multiplies the zero order equation of 
motion of the phase field (p. As a consequence, under the field redefinition 



we find 



0=0 + ^(f,Z>,C) 

So [0] = S [$\ - j d?x^gV 3 ) 5<f> 



(2.18) 



(2.19) 



In other words we are free to use the variable <fi instead of (p; however the first 
derivative correction with this choice of variable, Sl ven , differs from Sf" en by 



geven Qeven 

In other words the field redefinition fl2.18| ) induces the shifts 

A-A = o, / 2 -.f2 = o, h-h = H 



(2.20) 



(2.21) 



(where fi, and are the functions that appear in the expansion of Sf ven , 
see (|2.17|) ) For this reason, the dependence of all physical quantities - like the 
fluid constitutive relations - on f 3 is rather trivial, and easy to deduce on general 
grounds, as we will see below. 



2. While the fields a, \i and x are even under the action of time reversal, the 
fields and Q are odd under this operation. It follows that each of the three 
terms in (12.181) is odd under the action of time reversal. In other words the 



simultaneous requirement of parity and time reversal invariance simply sets W± = 
0. It follows that time reversal invariant superfluids have no non dissipative 
transport coefficients at first order. 



The corrections from ( 2.1C ) to the charge current and stress tensor Q2.14 ) in equi- 
librium are given by 
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In equations ( 2.22 ) and fl2.23| ) all the scalar functions fx, f%, fz and / have been treated 



as functions T, v and (C 69 ) 2 respectively. In obtaining ( |2.22|) we have used the zeroth 
order equation of motion for 0. 

= o 

to simplify the expressions presented above . 
2.5.2 Parity violating terms 

The most general parity odd contributions to the action are given by 12 

oodd Qodd I q 

*-> — &x ' '-'anom 

S°x dd = J ( gi & k CidjA k + T,g 2 e^ k Qd 3 a k ) + ^ j ada (2 2g) 
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13 where 



Cx is a constant and 



9i = gi{f,v,if>), g 2 = g 2 (f,u ) if>), 



A , C 2 

V = — , %) = — ■ 



(We emphasize that we have slightly changed notation compared to the previous sub- 
section. The independent variables for all functions in this subsection are T, v and ifj. 
The corresponding variables in the previous subsection were T, v and C 2 .) 

Note that ( J2.26 ) is automatically even under time reversal. The corrections induced 



by fl2.26| ) to the stress tensor and consistent charge current ([[27], see section 2.3 equation 



12 Our convention is | J XdY = j d 3 x y /g^e tjk X l djY k . 

13 The action for parity odd (non super) fluids ( |2 . 5| ) also contains the terms 

W=!f J Ad A + J Ada. 
But using the fact that d = Ai + di<j> and 



we can absorb Co in gi and C2 in 52 ■ 



2.16 of fl2]| ) in equilibrium are given by 
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where 
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Vl = e^ k F jk , Vj = ( : ' k f,,.. 



(2.29) 



The symbols for V$ and Vf have been introduced for notational convenience only; these 
vectors are determined in terms of the other quantities above by 
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-((C q ) l s,-v;) 



T 2 tp eq 

vi = ^—((c q ys 2 -vi). 



T 2 ip. 



(2.30) 



As we have emphasized, the formulas above determine the consistent current. The 
covariant current is obtained from the consistent current by an additional shift (see 



section 2.4 of [[R] for a review). We find that the one derivative contribution to the 
covariant current in equilibrium is given by 



Jo = -e CT {gi,pS\ + T g 2:0 S 2 ) 

t(/-eq\i yi (feq)i yi „ 1 

+ i(C q y(S igi ^ q +T S 2 g 2 ^ q ) 



(2.31) 



3. Constraints on parity even corrections to constitutive rela- 
tions at first order 

In this subsection we will determine parity even first order corrections to the superfluid 
constitutive relations both from the method of entropy increase as well as from the 
partition function of the previous section, and demonstrate their equality. 

Let us first consider the almost trivial case of parity even superfluids that also 
preserve time reversal invariance. As we have explained in the previous section, in this 
case W\ = 0. It follows immediately from this result that all non dissipative superfluid 
transport coefficients must vanish. Exactly this conclusion was reached in from 
the requirement of point wise positivity of the divergence of the entropy current in an 
arbitrary fluid flow. 

The study of time reversal non invariant superfluids is more involved. In this case 
the constraints from the local second law have not previously been analyzed. In this 
section we first present this analysis. We then study the constraints obtained from the 
analysis of the partition function. As mentioned above, we will find that these two 
methods yield identical constraints. 

3.1 Constraints from the local second law 

In this subsection (but nowhere else in this paper) we consider the non equilibrium flow 
of a superfluid on a (generically) non stationary spacetime. We continue to denote the 
background metric of our spacetime by G^ u . The background gauge field is denoted 
by A^. The variables of superfluid hydrodynamics are the temperature field T(x^), 
velocity field -u^x^) and the gradient of the phase field £ M = —d^cj) + A^. We often 
work in terms of the fluid dynamical field 



(C/)m = £m + A«V 



Note that, in equilibrium and at lowest order in the derivative expansion (C/)o — and 

(C/)i = 6 - A^i = Ci- 



We specify some additional notation that we will use extensively below. 




(3.1) 



Hi = T, H 2 = v, H 3 = (Cf) 2 



In words, P^ u projects onto the three dimensional subspace orthogonal to the nor- 
mal fluid, while projects onto the two dimensional subspace orthogonal to both 
the normal and superfluid velocities. a M and are the normal fluid acceleration and 
expansion respectively. V^ 1 is the 'Einstein combination' of the electric field and deriva- 
tive of the chemical potential that vanishes in equilibrium. Hi, H 2 and H% are new 
names for the three scalar hydrodynamical fields; note that H 2 is v = ^ rather than 
the chemical potential itself. Finally K is the term that is set to zero by the first order 
equation of motion of the Goldstone phase, while R is a combination of zero order 
thermodynamical fields that often appears in the formulas below. 

In order to analyze the constraints that follow from the local form of the second 
law, we follow the procedure described in section 3 of 0. Briefly, we first write down 
the most general onshell independent first order entropy current allowed by symmetry. 
We then compute the divergence of this current (this is mere algebra) and then use the 
equations of hydrodynamics, together with the corrected constitutive relations 



(here and refer to as yet unspecified one and higher derivative corrections to 
the constitutive relations) to re express this divergence as the sum of a linear form in 
onshell independent two derivative data and a quadratic form in onshell independent 
one derivative data. Point wise positivity of the divergence requires the linear form to 
vanish (this imposes several constraints on the entropy current). Once these conditions 
are imposed, the divergence of the entropy current is purely a quadratic form in one 
derivative data. We require this quadratic form to be positive definite. This require- 
ment further constrains the entropy current as well as the first order contributions to 
ti^ u and in a manner we now schematically describe. 

As we will see below, the quadratic form so obtained has the property that it 
vanishes when projected onto a subset of one derivative terms. In other words, all 



= (e + p)uV + pG"" + f^C + 7T 
J» = qu »-fe + f, 



(3.2) 



independent one derivative terms can be divided into y type 'entropically dissipative' 
terms and x type entropically nondissipative terms, and the quadratic form takes the 
schematic form 

Note that the structure of this quadratic form is preserved under x redefinitions 

but not under analogous redefinitions of y % . In other words there exists a well defined 
subspace of dissipative data but no definite subspace of nondissipative data. 

Positivity of the quadratic form described above requires that Aij is a positive 
matrix, and B im = for all i and m. The last set of constraints yield relations between 
otherwise apparently independent transport coefficients. 14 

In order to actually implement this process we need first to choose a basis for onshell 
independent data. As explained in || (see e.g. Table 3), at first order in the derivative 
expansion there exist 7 (4 dissipative and 3 non dissipative) onshell independent scalars 
, 7 (2 dissipative and 5 nondissipative) onshell independent vectors and 2 (1 dissipative 
and one nondissipative) independent tensors constructed out of thermodynamical fields 
and background fields. For the purposes of this section, we will find it useful to choose 
our onshell independent basis as follows. 

Basis of independent scalars: 

jjr^-, (u.dH a ), ((Cf).dH a ), a = {1,2, 3} 

The four of these scalars are dissipative (they vanish in equilibrium) while the remaining 
three are nondissipative (they are non vanishing in equilibrium, and do not cause 
entropy production). 

Basis of independent vectors: 



PT„, P^iCfhai, P£((fU va , P a »((f)„F» a , P» a d a H a , a = {1,2,3} 

The first two vectors are dissipative (they vanish in equilibrium) and the remaining 
five vectors are nondissipative. 

14 Assuming that the matrix A is positive definite, entropy is always produced whenever any of the 
y l are nonzero. It follows that all y l must always vanish in equilibrium. This observation motivates 
the following definition, utilized in Jl^] . Expressions that vanish in (arbitrary stationary) equilibrium 
are referred to as dissipative data. It follows from that entropically dissipative data is necessarily 
dissipative. However the converse is not necessarily true; it is possible for data to vanish in arbi- 
trary stationary equilibrium but yet be entropically nondissipative. We will see an example of this 
phenomenon later in this paper. 



Basis of independent tensors 



(J = P a P 13 



T (Cf) 

'ill/ 



pa pP 

(I V 



V (C/)/J + Vf>{Cf)a ~ P X ^x((fh)V a p 



tive. 



The first is dissipative (it vanishes in equilibrium) while the second is nondissipa- 
In this paper we wish to constrain the equations of superfluid hydrodynamics 



presented in a 'fluid frame' (see [28 for an explanation of what this means). Throughout 
this paper we will further restrict our attention to fluid frames with /^ ss = (again see 
for definitions). This choice still permits the freedom of field redefinitions of the 



temperature and normal velocity fields (as well as field redefinitions of the superfluid 
phase, as we will exploit later in this paper). Even though we work specifically frames 
in which Hdiss — our final results may easily be lifted to an arbitrary fidiss 7^ frame 
using the frame invariant formalism of 0. 

The most general form of the entropy current, consistent with the absence of linear 
two derivative terms in its divergence was determined in || (see equation 3.19 ) and 
takes the form 



7M _ , J(J. 



u v ^ v 
T 



a 

where QT = f(u»(QY ~ u"{CfY) 



(3.3) 



The divergence of J^ an was worked out in [28|, || (see for example, equation 3.9 
and recall we work with Hdiss — 0)- 



V J* 1 

v (i" can 



^'T ll [^)+fV, + (u,f)(u.du) 



(3.4) 



The RHS of ( |3.4|) is given schematically by 
(one derivative correction to constitutive relation) x (entropicallydissiaptive data), 



15 We will now rewrite the RHS of ( |3.4| ) as a quadratic form in the basis of independent 
dissipative one derivative data chosen above. In order to achieve this we need to rewrite 



15 Note that the one derivative expressions that appear here are always entropically dissipative, as 
contributions to changes in the proportional to these one derivative expressions yield quadratic terms 
in entropy production. 



all the y type terms in ( |3.4| ) in terms of the independent basis of dissipative scalars, 
vectors and tensors listed above. To achieve this we use the equations of motion 



(QUtfUZ . 5 = _ T n - „mYM _J^L_ MXCf) 2 (3 ' 5) 

(C/) 2 3 " 1 ^ (C/) 2 e + P 2(C/) 2 



we find 



V 7 M 

v A 4 " can 



+ 



T 



[u.dT) + (j.u)(u.du) + 
if 



PT 



mc,f) 

V (C/) 2 



(j-(C/))(K(C/)) 
(C/) 2 

- h 



i 

T 



nc f y 



e + P 
(C/) 2 3 



J^2 



[IV 



[u.dT) + (J.u)(u.du) - 



TT^P, 



[iv 



2T 



e 



+ I TV^P^ - -TT^P, 



/ ^•d(C/) 2 

V2T(C/) 2 



+ 
+ 



V-(Cf) 
(C/) 2 

K 
e + P 



(j.(C/)) + P(« U (C/),^) + (1 - //P) ( TT^ - -TT^P^ 



-U 



71 



[IV 



T 



+ V TT^P 



-7T 



P 



[IV 



( (fU a »Pv^(( f ) 

ncf) 2 



+ (j» + Ru a 7r a »)P^V--a 



' [IV 



= g G a (u.d)H a + 6 4 (YM) - mi {^ V T{Qf" ) + ~ \°^ V 

(3.6) 
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(3.7) 



The last line of ( |3.6|) is the final result of this manipulation. It expresses the divergence 
of the entropy current as a linear sum over the four dissipative onshell scalars and two 
dissipative onshell vectors and one dissipative tensor listed earlier in this subsection. 
These expressions appear multiplied by frame invariant linear combinations of 7r Ml/ and 

The frame invariant quantities & a and QJ a will be used extensively below. For later 
use we will find it useful to regard these quantities as linear functions of tc^ and 
i.e. 

e a = e a (TT^,f), W a = X a (n^,f) (3.8) 
The divergence of the 'new' part of the entropy current, J£ ew (see ( |3.3| )) is given 

by 

V P 

v A 4 " new 
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(C/),V M Q^ 



(3.9) 



.dH a ) (u u V u Q^) + ((Q).dH a ) 



where Q^ u was defined in (|3.3| ). 

Using equations of motion we can express (u^VuQ^), ( (c f ) 2 — ) an d P^(V a (5 aiy ) 
in terms of the onshell independent basis scalars of this subsection ( spanned by 



i u 



.dHa), ((Q).dH a ), 



((f) 2 



and vectors (spanned by P^K, P M °(C/)/3^, P M °9 Q P a ). 
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(3.10) 



From equations (|3.6|) , ( |3.7| ), ( |3.9|) and ( |3.10| ) we conclude that, no matter what form the 
fluid constitutive relations take, the divergence of the entropy current cannot contain 
any expressions of the form ((Q).dH a ) 2 or (P^d^HaduHb)- In other words the scalars 
(Q).dH a and the vectors (P^d^Ha) are nondissipative. It follows that the positivity of 
(V^Js) requires that the divergence contain no term linear in ((Q).dH a ) or (P^d^Ha) 
(see e.g. P for repeated use of similar arguments.) To ensure this tt^ and j M have to 
satisfy the following conditions. 
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where M is a 4 x 4 matrix of dissipative transport coefficients in the scalar sector and 
TV is a 2 x 2 matrix of dissipative transport coefficients in the vector sector. 

Equations (|3.11| ),( j3.12j ) are the main result of this subsection. It expresses the 



equality type constraints that follow from the local second law. Once fl3.11|) , (|3.12| ) 
are satisfied the final expressions for the divergence of the entropy current takes the 
following form. 



- E c ° + s ( 1 ~ 7^) ^ + / ^ + (u.dH a )(u.dH b 

- f(V.(( f ))J2ca(u.dH a ) 

a 
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+ f 

(3.13) 

The positivity of this quadratic form imposes additional inequality type constraints on 
transport coefficients that we will not further explore here. 

3.2 Constraints from the partition function 



In this subsection we now reproduce the conditions ([TTI]),([JT2]) using considerations 



independent of those of the previous subsection. The procedure we adopt is very 



similar to that described in ||12|| , and we describe it only briefly, highlighting only those 
elements of the analysis that are unique to the superfluid. 

The starting point of our analysis is the expressions ( [2.22 ) and ( p.23| ) which rep- 



resent the first order for the corrections to the stress tensor and charge current that 



follow by varying the local action for the Goldstone mode w.r.t. the metric and back- 
ground gauge field. Once we substitute in the solution for the field ^(x), according to 
its equations of motion, ( 2.22| ) and fl2.23| ) yield first order corrections bT^ v and 5J M to 



the values of the stress tensor and charge current in thermal equilibrium. 

From the hydrodynamical point of view, 5T^ V and 5J^ are the first order contri- 
butions in (13.21) once we substitute 



T(x) = f(x) + T 1 (x), fi{x) = p.{x) u^{x) = m m + (3.14) 

into those expressions. Here Ti(x), /ii(x) and Ui(x) a are the first derivative corrections 
to the equilibrium configurations of temperature, chemical potential and velocity. 

Upon substituting ( |3.14j) into ( |3.2|) we get first derivative contributions of two 



sorts. First we have the corrections to constitutive relations evaluated on the zero 
order equilibrium configurations W ,v (f,fi,u^ 1 , {(, eq Y) and j M (T, {C, eq Y) ■ Second 
we have contributions from terms proportional to Ti, /ii and Ui when ( |3.14| ) is plugged 
into the perfect fluid constitutive relations. Contributions of the second sort, however, 
precisely cancel out in the frame invariant linear combinations & a (a = 1 ... 4) and QJ a 
(a = 1 ... 2). In other words 

(3.15) 



7r MV and j M on the RHS of ( |3.15|) are evaluated on the zero order equilibrium configura- 



tions). In the general formulation of hydrodynamics, however, it is precisely the frame 
invariants that appear on the RHS of ( |3.15| ) that are expanded in the most general 
symmetry allowed constitutive relations (see e.g. ) 

e a (n^,f) = a am S m (a = 1...4), (m = 1 . . . 7) 
W a (n^,j»)= lam V£ (a=1...2), (m = 1...5) 

where S m and VJ^ are the independent one derivative scalars and vectors and the 
coefficients a am and 7 am are arbitrary functions of the scalars T, /i and 

a am and 7 am are the constitutive coefficients we wish to constrain, and this is 
achieved as follows. In the LHS of ( |3.15| ) we substitute the expressions (|2.22|) and 



( |2.23D ) for ST^ and 5J M . This determines the LHS of ( |3.15| ) completely in terms of 



the functions fx, ji and that appear in the partition function. In the RHS of ( |3.15| ) 
we substitute ( |3.16|) , and evaluate these expressions in equilibrium 

T = f, n = ( = C q - 



Under the last substitution those of S m and V m that are dissipative vanish. The 
non dissipative one derivative scalars and vectors evaluate to geometric expressions. 



Equating the coefficients of these expressions we determine a am and j am for those 
values of m that correspond to non dissipative terms. In other words this procedure 
completely determines all non dissipative transport coefficients. 16 

In the rest of this subsection we implement the procedure described above to ex- 
plicitly determine all nondissipative transport coefficients in terms of the three free 
functions fx, fi and f% that enter the local action for the Goldstone field. We demon- 
strate that our results agree exactly with ( |3.11|) , (|3TT2| ), obtained from the local form 



of the second law, once we identify the three unknown functions ci, c 2 and C3 in the 
entropy current of the previous subsection in terms of the functions in the partition 
function according to 

A L?k A L?k L?A 
Cl fT + TgT , C2 fT + Tdy i c 3 Td( - 2 

We will also demonstrate that the identification fl3.17| ) may be argued for directly 
by comparing the thermodynamical entropy in equilibrium with the integral of the 
equilibrium entropy current over a spatial slice. 

It will be useful in the computation below to note that P^ v and P^ v are given by 

ij ~ 9iji "ij — 9ij ~ J7eq\2 

We turn now to the explicit computation, starting with the vectors. 

QVWV,^) = 5J 20 (5T^,5J M ) = 

SJuiST^SJJ = P i3 (Sj j + Ru°5T^j 

= Pij (5P - Re- a A 5J j ) 
= (1 - fifyPijSJi 

= Py^*(l - jlR) (kd k f + kd k V + Idkfi 

= - P^g* ( f -ld k f + hd k v + ld k f 3 ) 
\T T T J 



(3.18) 



16 There is an important subtlety here. All of the operations described above may only be performed 
in equilibrium, i.e. once we have solved for (C 51 ?)^ as a function of background fields and substituted 
this back into the partition function. We implement our programme without explicitly solving, simply 
by treating (C eq ) fl (x) as formally independent of the other background fields, except for those local 
combinations of ^ e i)^ that appear in terms of its equation of motion and derivatives there off. The 
reason for this that the expressions for £ M as a function of background fields is highly nonlocal. 
The only situation in which cancellations are possible between local expressions in (£ e| J)^ and local 
expressions in the background fields is when we get derivatives combining with (£ ei ?)^ in the form of 
the 4> equations of motion. 



The last line of ( |3.18| ) exactly matches ( |3.11| ),( |3.12| )) upon using the identification of 
the parameters (|3.17| ). 

We turn next to the scalars; let us start with ©4. 



(3.19) 



6 4 (£ZV, s j M ) = (j.c) + R{u^ v ) + (1 - nR) [ ^ v p»v - 2^ vp ^ 



- (1 - jlR) [A(C q .dT) + f 2 (C q -dO) + f(C q -dh) 
T(l-fiR)J2fc b (C q -dH b ) 



In the last step we have used ( |3.17| ), and have obtained manifest agreement with 

Next we shall calculate the remaining three scalars & a , a = {1,2,3}. The 
algebraic manipulations here are a little more involved than in previous cases, and we 
provide some details. 
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The first line in Q3.20D can be evaluated as 
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In the last step we have used ( |3.17| ). 

The second line of ( |3.20|) may be evaluated as follows 



1 ds 
~s~dH a 

1 ds 
~sl)H a 
1 ds 
^dH n 



7T 



\iv p 



2T 

^{C q .d)f + ^{C q .d)v + L {c ^ d)h 



J2fc b (C q -dH b ) 



(3.20) 



(3.21) 



(3.22) 



In the last step we have used ( |3.17| ). 

Finally we evaluate the last three terms of ( p.20|) together. 
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(3.23) 



In the last step we have used (|3.17| ). 

Combining (|3.21|) , (|3.22|) and ( |3.23|) it is straightforward to verify that the expres- 
sions for & a , a = {1,2,3,4} as derived from partition function in this subsection, 
match exactly with (|3.11|) , (|3.12|) . Note that both methods leave dissipative contribu- 
tions to constitutive relations completely unconstrained. 

3.3 Entropy from the partition function 

In this subsection we will explain how the nondissipative part of the entropy current 
of the superfluid may be read off in a rather direct way from the partition function. 
Our analysis is largely structural, and applies equally well to normal (non super) fluids. 
However our presentation applies only at first order in the derivative expansion. 

For any system the entropy St in equilibrium may be evaluated from the logarithm 
of partition function W = In Z via the thermodynamical relation 
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We will now rewrite this expression in terms of the goldstone action that generates the 
partition function. Let this action take the form 



S = J ^Cd 6 x 

and also suppose 

Now we can think of the partition function as 



W=S(T,fi,T a i ,C q ) 



Using the simple rescaling of the time coordinate employed in subsection 2.3.1 of [[12 
one may show that 
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so that 



S T = W + 



dW 
dT 



y, V-G 



d s y 



To 



WT® 
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(3.25) 



(3.26) 
(3.27) 



This expression may be expanded to first order in derivatives employing 



C eq = ^ + Cf 
T 

To = ( T o)perf + 
J° = J° P erf + SJ° 



(3.28) 



where, from (|2.14j ) 



CO 



perf 



e-fe-^AZ-fAoaX 



is-eq 



Jperf = e- a q - f(e-^A + a*CF 

5T ° is defined in ( ggg ) and 5J° is defined in (ggg ). 
Using the Gibbs Duham relation 

P + e - qfi 



(3.29) 



T 



we find that 



d 3 yy/gs 



+ <Ty 



(3.30) 



(all proportional to / cancel out at zero order in the derivative expansion). 
Now let us recall that 

ST^ = (7r^) + (Tf?)p er f 

where (7r^) refers to 7r^ evaluated on the zero order equilibrium solution and {T^Yperf 
refers to the one derivative correction in Tj^ from the first order correction to the 
equilibrium solution. Similarly 



It follows that 



S T = J d 3 yV^G 
+ / d 3 y 



^^(T^perf 



Tn 



su — 



TC\ q - 



T 
f 



perf 



- Mf 



(3.31) 



so that 



Sn 



d 3 yV^Gsu° 



+ d% 



To 



fCf - ^> - 05 f 
T 



(3.32) 



where su° in ( |3.32| ) refers to the entropy evaluated on the first order corrected solution. 
In going from (|3.31|) to (|3.32|) we have used the fact that the frame invariance (see 
for a definition and extensive discussion of frame invariance) of the canonical entropy 
current 



"can 



U,,TT' 



implies that 



sir 



iperf 



It follows from (13.321) that 



S T = j d 3 yV^G 



+ 



T 



perf 



T 



1 



0. 



T 1 



Comparing with 



we conclude that 



J s 



can ' " new 



J d 3 yV^GJ° new = f dZyV^G^Cl 9 



(3.33) 



(3.34) 



(3.35) 



In other words the integral of J® ew matches with the first order correction to the 
Goldstone action. ( |3.35|) is the principal formal result of this subsection. It expresses 
a very simple relationship between the correction to the canonical entropy current of 
our system and the first order correction to the partition function. 

To what extent does ( p. 35 ) determine J£ ew 1 The most general first order correction 
to the entropy current takes the form 



"new 



(3.36) 



where S u and S$ are first order scalars while Vf is a first order vector. Notice that, 
to first order, X M = + Vf 1 is orthogonal to u. It follows immediately from this 
observation that 

X° = -a,X i 



Plugging this relation into ( |3.35| ) we conclude that the contribution from X M to the 
total entropy is not Kaluza Klein gauge invariant and so must vanish (see |T2| for a 
discussion on related issues). It follows that and Vf vanish in equilibrium. Upto 
dissipative corrections, therefore, it follows that 



°new ~ ^u" 

Now comparing with (|3.35|) it follows that 



(3.37) 



J d 3 y^(S u -C e 1 q ) = 



so that 

S u = C eq + total derivatives (3.38) 

Let us now turn to the case at hand. C eq was listed in ( [2.17] ). It is easily verified 
that there exist no total derivative scalars at one derivative order. Consequently we 
conclude that 

S u = j(C-d)T + j{C-d)p - / 3 V, (IA + dissipative 



It is not difficult to verify that this expression, together with ( |3.37|) , agree exactly with 
( p.3| ) in equilibrium once we employ the identification of parameters ( |3.17| ). 

In summary, the positive divergence entropy current - which we determined earlier 
in this paper - is also uniquely determined by comparison with the partition function 
for parity even superfluids at first order in the derivative expansion. 



3.4 Consistency with field redefinitions 

We will now verify that the dependence of the constitutive relations and entropy current 
of the superfluid on f 3 is consistent with the transformation ( |2.21| ) of / 3 under the field 
redefinition ( [2.18| ). 

Recall that the stress tensor and currents of our system take the form 



fJbV 



T» v = (e + P)uV + PC" + feC + 7T 
J» = qu »-fe + f (3.39) 

J a ~ ^ can J new ~ SU ^ Tj, ^j^ + J new 



Substituting the field redefinition (|2.18|) into this equation and setting 



S^x 1 ) = h{x l ) 

(recall h is a function only of space) we recover a new form of the stress tensor and 
currents 



T» u = (g + P)u^u u + PGT + feC + r 

r = qu» - f~e + ? (3.40) 

u u n> MU 

Js ~ Jean ^ new ~ SU ^ rn V j ^new 



with 



7T 



71 



f- 



d(e + P) dP df 

OX ox ox 



-2£.V (4) /i) 



dx dx 



-2£.V (4) /i) + fG» a V^h 



(3.41) 



new 



T 



-v(f-f) 



All Greek indices in (|3.41|) and ( [3.40 ) run from 1 ... 4 and are raised and lowered with 
the full four dimensional metric G^ u . x derivatives in Q3.41 ) are taken at fixed T and 
v. In deriving last equality in ( |3.41| ) we have used the first law of thermodynamics. 

/ 

de = Tds + vdq — — dx 

We will now independently verify that our final answers for J£ ew and the constitu- 
tive relations have this symmetry. To start with recall that, from ( |2.21| ) and (|3.17 ), 



1 dh 
TdW a 



(3.42) 



It follows immediately from (|3.42 ) that the expression for J£ e 



(see (|3.3| )) transforms as predicted by the last of Q3.41] ). 

We now turn to the verification that our results for transport coefficients, ( |3.11| ),( |3".12| ), 
transform as predicted by ( 3.41 ). The algebra involved in a direct verification is 
formidable, so we will content ourselves with an indirect check. We first recall that we 
have already verified (see ( |3.15| )) that 



& a (5T tl „,5J ll ) = e 
W a {5T tlv ,8J lM ) = W 



(3.43) 



in fact this equation formed the basis of one of our two methods of determining consti- 
tutive relations. It follows that if we can show that 5T flu and 5J^ obey 03.41 ), then the 
same will be true of (|3.11|) , (|3.12|) . (Recall 5T^ U was the first order shift in the stress 
tensor arising from first order corrections to the Goldstone action; 5J^ was similarly 
defined.) We will now check that this is indeed the case. In order to do this we first 



simplify the ( |3.41| ) specializing to the case of stationary equilibrium 





(j -J ) = -2e° 



dcf q + fif) 



(C eq .d)h 



(C q .d)h 

(f -f) = -fV*h-2(C*-dh)^ i (C 







dv 


m 



and 



{n 00 - 7r 00 ) = -2e 



2a 



2e 
2e 



2(7 



da 



d (BP „ 
T— - P 



80 V ~ dT 



{C q -d)h 

(C q -d)h 



2dT 2 



2 2a 



-T z e 



\- 




dT 


8) 



(C q .d)h 



(3.44) 



K-^) = -A (/-/) 



(yi - = 2{C, eq .dh) 



d0 



+ f[(( eq ) i V j h + (C 9 ) j V i h] 

(3.45) 

Where each of the scalar thermodynamic functions are evaluated on the zeroth order 
equilibrium solution 

t = f, u = d, (C/), = cr 

In obtaining (|3.44j) and (ft.45|) 

'e + P + fi 2 f 



dP 



) dT + T{q + uf)dv- f -dQ 



In those equations all spatial indices are raised and lowered by use of the spatial metric 
Qij (all the free indices will run from 1 to 3). 

We now turn to the explicit expressions for ST^ and 5J^ listed in (|2.22|) . Substi- 
tuting 

~h = h + h 



(see ( |2.21| )) in those expressions we obtain immediate agreement with ( |3.44| ) and ( 3.45|) . 
This completes our verification. 



4. Constraints on parity violating constitutive relations at first 
order 



In this section we use the partition function to derive constraints on parity violating 
contributions to constitutive relations by comparison with the local goldstone action 
( |2.26D . As in the previous subsection, we find perfect agreement with the constraints 



obtained from the local form of the second law. It turns out in this case that the 
second law analysis has already been performed, in full detail, in |J. We begin this 
section by reviewing the results of ||, before turning to a re derivation of those results 
by comparison with (|2.26 ). 



4.1 Review of constraints from the second law 
4.1.1 Basis of Frame Invariants 

As we have seen above, the constitutive relations are an expansion of frame invariant 
combinations of 7r MV and j M in terms of independent one derivative scalars, vectors and 
tensors. Before even specifying the constitutive relations, we must first specify a basis 
of frame invariant expressions that we will expand in this manner. In the previous 
section we choose to work with the frame invariant scalars & a and frame invariant 
vectors QJ a . A different choice for frame invariants was made in ||; in order to ease 
comparison with the results of that paper, we will adapt that choice in this section. In 
this subsubsection we describe the basis of frame invariants used in ||. 
Let 

Si = vr^/V C/ 2 S2 = C/-vr-C/ (4.1) 

S 3 = U-7T-U S 4 = U-7T-(f 

S5=U-j S 6 = C/-J 
s 7 = ~^diss 

u^ a P\ ^ = (C/)x°n 



V 



V3 = P U af 



2 



Throughout this paper s 7 = —Hdiss = ( Hdi SS was defined in |28|1 ). However we will 
retain in all our formulas, in order to permit easy adaptation of our final results to 
frames in which Hdiss 7^ 0. 

(C/r 



Following || we define the row vectors 

S = ^Si S2 S3 S4 S5 S6 S7 
v = (v X v 2 v 3 



We also define the matrices 
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\-i / 
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p + P T 
and the Aj's i?;'s, Cj's and iYj's defined as follows. 



P^d^v 
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d -[log(s)], B 3 = -^[log(s)} 



P 



dv 
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(4.3) 



(4.4) 



(4.5) 



In terms of (IO 



the frame invariant scalar, vector and tensor combinations 



of tt^ u , and Hdiss are given by the row vectors 



sA s 



(4.6) 



By scalars, vectors and tensors we mean expressions which transform as spin 0, ±1 and 
±2 representations of the SO(2) symmetry that is left invariant by the two vectors 
and £^ at each point in spacetime. 



vector 



definition 



(f -vft)) 



T^V, (£) 



dual parity odd vector 
evaluated in equilibrium 







-TV? 



I P cr V i 
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Table 1: Independent fluid vector data. Here for m=l,2,3,4,5 are independent vectors 
in equilibrium defined in (|2.2S| ) 



pseudo scalars 


definition 


In equilibrium 


<Si 






s 2 




"1 + £o<S2 


pseudo tensors 


definition 


In equilibrium 


fr 







'2 




won't need 



Table 2: Independent fluid scalar and tensor data. Here S m for m=l,2 are independent 
vectors in equilibrium defined in ( 2.2S| ). a u and are the shear tensors for normal and 
superfmid velocity respectively and *Ouv = e flpal3 n p f a i- 



4.1.2 Constitutive Relations 

We have 4 frame invariant scalars, 2 frame invariant vectors and one frame invariant 
tensor. The most general symmetry allowed parity odd first derivative constitutive 
relations take the form 

2 



i=l \i=3 
2 4 



(4.7) 



,i=l 3=1 



with T, T, V, V, iS and S a basis of onshell independent SO (2) invariant tensors, 
vectors and scalars given in table 



4.1.3 Constraints on constitutive relations from the local second law 



Notice that both pseudo tensors that appear in ( |4.7[ ) are nondissipative. Further, the 
five pseudo vectors Vi % = 3 ... 7, are nondissipative. It may therefore come as no 
surprise to the reader that f| was able to use the principle of local entropy increase 
to determine ki m (i = 3 ... 7 and m = 1 ... 2), together with (i = 1 ... 2 and 
j — 1 ... 4) in terms of two free functions that appeared in the parameterization of the 
entropy current. These two functions were called a 8 and a w in ||. The results of f|] 
were presented in terms of a 8 and cxio and four additional auxiliary fields which were 
determined in terms of a 8 and o"io by the relations 17 

(J 3 = -T-^(a 10 - va 8 ) 

d 

a 4 = cr 8 + Cv + 2h - — (cr 10 - vo 8 ) 

g d " (4.8) 

^5 = -^—(^10 - va 8 ) 

2 

a 9 = 2v{a w — ua 8 ) CV 3 — 2hv 2 + s 9 

3 

In terms of all these fields, it was demonstrated in || that point wise positivity of 
the the divergence of the entropy current determines 

V = 0, k m2 = (4.9) 



and 



K31 = —RTcr 3 — Tdrcr 8 
K41 = —RT 2 a 4 — Td v a 8 

k 51 = -RT 2 a 5 -Td, 4; a 8 (4.10) 
« 61 = -2RT 3 a 9 + 2TVi 
«n = -RT 2 a w + 2Ta 8 + CTz/ + 2/iT 

2HZV9 _ 2*710 _ 2cJ 3 _ 2T 2 K 3 a 9 -2T0-4 - 2T 2 K 2 cr 9 -2Ta 5 - 2K l T 2 a 9 



lJ \- Q T^ 1 ~ + d T V8-K 3 Ta w d u a 8 - K 2 Ta w d^as - K x Tg 1q 

(4.11) 



17 All terms in (4.8) proportional to the constant h were omitted in ||. The reason for this is that 
H assumed that the entropy current was gauge invariant. As explained in |]l2| this does not seem 
to be physically necessary as long as the divergence of the entropy current is gauge invariant. This 
allows the addition of the new term proportional to h in (121), which allows for a slight modification 
of the results of ||, captured by the shifts described below. As we will see later, the requirement of 
CPT invariance forces h to vanish. 



4.2 Constraints on constitutive relations from the Goldstone action 



As in the previous section, we use the Goldstone action to constrain transport coeffi- 
cients as follows. All constraints follow from the analogue of ( |3.15|) 

vfA^ST^SJ^ = vfA^Tr^,^) (4.12) 

The LHS in this equation may be determined in terms of the functions g\ and g 2 in the 
Goldstone action using ( |2.28| ). The RHS of the same equation is simplified using ( |4.7| ) 



under the substitution T — > T, \i — > /t, Q — > Q eq . Under the last substitution, the parity 
odd first derivative vectors and scalars evaluate to geometric expressions. Substituting 
these relations into the RHS of ( |4.12j ) and equation coefficients of independent vectors 
and tensors yields an expression for all non dissipative transport coefficients in terms 
of the functions g\ and g 2 . Using Eq.( |2.27| ) one obtains 



t 

f 2 ,^ „ f 2 



= T 3 (-ud fgi + d f g 2 )V{ + ^{vdogx - d fi g 2 )\q - j^( 9 2 - 2g x i>)V^ 
V3 = Kf 

= f(fd f9l VI - ±d*9i Vi - J^i^i Vi + 9 2 T Vi + (C q )% eq9l Vi)) 
Si = Tr^iv = 

s 2 = tt^C/^-C/ = -l{C q ) 2 {d^ q9l 8 X + T d feq g 2 S 2 ) 

s 3 = u-tt-u = f(fd f g 1 - 2ijj eq d i>eq gi)S 1 + f T (fd f g 2 - 2^ eq d^ q g 2 )S 2 

s 4 = u-7r-Cf = {f 2 (g 2 - 2 9l ) - 2{C q fvd^ q9l )S 1 - S 2 T D(g 2 f 2 + 2{C q fd^ q g 2 ) 

+ 2Qe- ff 5 2 T 3 - C-Ale~° (A S 2 + 3S 1 ) 
o 

s 5 = u-3 = -{dpgxSt + T d g 2 S 2 ) 

(C eq ) 2 - (C eq ) 2 1 

s 6 = Cfj = 2T(g x + ^>d^ q g 1 )S 1 + TT (g 2 + ^>d^ q g 2 )S 2 + -CA e^ (A S 2 + 2Si) 

1 z 1 z * 

S7 = —fJ'diss = 

(4.13) 

Now using Eq. ( |4.7|) one can find out the transport coefficients in terms of partition 



function coefficients g x , g 2 as follows 

fj = 0, k i2 = /or z G (3 to 7) 

f ((-i>aT + e + P)<9 t 0l + qfd f g 2 



«31 



«41 



«51 



P 



e 



f ((-%f + e + P)d gi + gfa^ 2 
PT^ 

f ((-%T + e + P)d^ q9l + qfd^g^ 



P + e 
2T 2 / 

«6i = [-g2{-2vqT + e + P) + 2g x v{-vqT + e + P) 

Cu 2 f 2 (3p-2Pqf + 3e) 4dgT 3 



+ 



3(P + e) P + e 

P / * A * \ Cvf(2p-vqf + 2t) 



«n = — p (<MT + 2 gi (-uqf + e + P)) + 



2(P + e) 

(4.14) 



Similarly, the transport coefficients /3y in terms of partition function coefficients gi, g 2 
as follows 

4Pfz>(-g 2 + qi z>) 2(-g 2 + 2 gi u) r , 2 f 2 (-3P + 20qf-3e) 4gT 3 
— Pll ~~ i / r O ; — r h La- 



^ ^ 3(C e «) 2 (P + e) ^C^^P + e) 

_ P(-g 2 + 2gj£) of (2P - vqf + 2e) 

Pl2 ~ n eq + ^ 2(c^) 2 (p+6) 

-/3 21 = -2T{-vd f g x + d t g 2 ) - 4z>T 2 K 3 (-g 2 + g x v) - ^CK 3 f 2 v 3 - AC x K 3 f 2 
-fa = d f9l - K 3 f(-g 2 + 2g x v) - \fK 3 fv 2 

-fa = -2T{-vd g x + d g 2 ) - 4uf 2 K 2 (-g 2 + 9l u) - ^CK 2 f 2 tf - AC\K 2 f 2 
-fa = dp 9l - K 2 f{-g 2 + 2g x v) - l -CK 2 fv 2 

-fa = -2T{-vd^ eq g x + d^ q g 2 ) - AK 1 0T 2 (-g 2 + g x v) - ^CK x f 2 u 3 - AC x K x f 2 , 
-fa = d^ q9l - K{f{-g 2 + 2 9l u) - l -CK x fu 2 . 



(4.15) 



In equations ( |4.13| ), ( |4.14 ) and (|4.15|) the functions g x , g 2 and all the other thermody- 
namics functions (like e, P, q etc) as arbitrary functions of T, v and ip eq . 
If we make the substitution 

01 = o- 8 + h, g 2 = -o 10 + 2z>o 8 + hjv 2 + 2hv. (4.16) 



and introduce the auxiliary fields 03, 04, 05 and <7g which are written in terms of erg and 
cio in ( [4.8|) then our results for nondissipative transport coefficients agree precisely 18 

with o, flog), (grip. 



4.3 Entropy 

As in the parity even case, we may determine the parity odd contribution to the entropy 
current by a simple direct comparison with the the partition function. The relevant 
equation here is 

W° dd + W anom = = J d 3 yV^G [P{5,P consistent - 6J° covartant ) + J° Sn J 

= J rf»yV=G [u8J° shift + J° Sn J (4.17) 



The term in (^4.17|) proportional to 5J^ hi j t has its origin in the fact that ( |3.27|) is correct 
when J° is taken to be the consistent U{1) current. On the other hand the canonical 
entropy current of hydrodynamics is defined in terms of the covariant U(l) current. As 



explained in |12| these two currents differ by the shift 



fsun = ^ vpa A v 7 p(J . (4.18) 



The contribution of this shift to the RHS of ( |4.17| ) evaluates to 



d 3 yV^GC>5J° hlft 



C 
~6 



- J d 3 y^ue tjk (^A i d j A k + ^A A i d j a k + ^A a i d j A k + Alaidja^ 

% I <Py,/gi>j jk (A i d j A k + -4,^0^* + -^—A (a.(C 9 )Si - a.V 3 ) 
A J V ^ 2T 2 i> eq 



+ -^L(a.(C«)S 2 -a.n) 



(4.19) 



18 We also need to make identification sg = 2Ci, as will be clear below. 



Taking this contribution to the LHS of eq.( [4.17| ) we find 



j new 



J d 3 yV^GJ; 
1 dd + W anom - J d?yV^G5J : 



o 

shift 



I d 3 yy/g\ 0i Si + T g 2 S 2 + 
J \ T 2 y 



(CX-^A 2 )(a.(C)S 2 -a.V 4 ) 
' CvA (a.(C q )Si-a.V 3 ) 



2T^ eq 



(4.20) 



In rest of the section we will use (|4.20Q to constrain the new part of the entropy 



current. The most general form of the first order entropy current is given by 



+ aiV^ + a 2 V 2 CM + q [a 3 {u-0 + <h(B-Q] + he^ Xr7 A u d x A a 

where h is a constant (4.21) 



Since the first term proportional to o\ is a total derivative, it is not determined. The 
term proportional to ol\ and a 2 is also undetermined as and V 2 ^ both are zero at 
equilibrium. We now evaluate Q4.21 ) in equilibrium. Using Table 1 and Table 2 and 



u° = ^((a.(C«))S 2 -(a.V 4 ))., 



V c/ = 0, V c0 ' T = -a i V ci ' 1 where I G (1 to 7) 
e a 

B° = -j^yi ((a.(r))5i - (a.V 3 )) - j^y 2 {(a.(0)S 2 - (a.V,)) 

\e^^ v F pa = <*{S X + A S 2 ) + ^(( a .(C*))5x - (a.V 3 )) + ^L(( a .(C*))& - (a.V 4 )) 
+ e a {a.V 2 ) 

e°» Xr7 A u d x A a = e~°e ljk [AdjA, + 2T va i d i A k + T^a^a* + d { (T Q v aj A k )] (4.22) 



where Vi and Si are listed in Eq. [2.28| . Now using the fact that (( eq )i = Ai + 

e ilk [AdjA k + 2T Oa i d j A k + T^a^a^ 
^ k [(C%d,{C q )k + 2T ua l d J (C q )k + Tya^ak] 

= J V9{Si + ^2T z>((a.(r))5 1 - a.V 3 ) + ^ T y((a.(C))S 2 - a.V 4 )) 

(4.23) 

we obtain 

J d 3 xV^GJ° Sncv = J rf 3 x v ^(e <T Ta 3 (a.y 1 ) + (a 8 - a 4 )(a.V 2 ) - fe°a 5 (a.V 5 ) 

+ as^ + A S 2 ) + hSx - (a.(C e<? ))( - \^e° S 2 + a 4 (5' 1 + A S 2 )) 

+ J^y(~ f eVl ° + asA ° + c ^T Q v){{a.{C q ))S 1 - (a.V 3 )) 
1 p 2ct . 

— TV 9 - Te** 10 Ao + a s A 2 + hT 2 2 ) {(a.(( e «))S 2 - (a.V 4 )) 



(C e 9)2 v 2 



(4.24) 



It is convenient to introduce the following redefinitions 

a 3 = -fd f X, a 8 -(T 4 = d X + Y, a 5 = -d^X + Z. (4.25) 

Now using 

d k X = d f X d k f + d,X d k u + d^X d k ^ eq , (4.26) 
the first line of the Eq. f4.24j can be rewritten as 



J tfx^g^Ta^a.VJ + (<r 8 - a,)(a.V 2 ) - fe*a 5 (aV s f 
J <Pxy/g(T P k a i (( e *) j d k X + Y(a.V 2 ) - fe"Z{a.V 5 j) 

J d 3 x^( - ToXe^iO^a, + T Xe« + Y(a.V 2 ) - Te"Z(a.V 5 f 
J d 3 ^( - T XS 2 + T X-^((a.(C q ))Si - (a.V 3 )) + Y{a.V 2 ) - fe°Z(a.V 5 )^ . 



(4.27) 



So we obtain 

J d 3 xV^GJ° Snevr = J d 3 x^( - T XS 2 + a 8 (S 1 + A S 2 ) + fa 

+ (T X - fe°a l0 + a 8 A + 2hT is)-^((a.(C q ))Si - (a.V 3 )) 

+ J^fC^T V 9 - fe°a 10 A + a s A 2 + hT^) {(a.(C q ))S 2 ~ (a.V 4 )) 

- (a-(C e, ))( " ^ 3 e CT S 2 + 04(^1 + 4>S 2 )) + F(a.y 2 ) - T e ff Z(a.y 5 ; 

(4.28) 

Now using (|4.20|) we obtain 
Y = Z = 0, a 3 = a 4 = 

0-3 = -T% (<t 1q - z>a 8 ), 0-4 = a 8 - do(<J W - z>o- 8 ) + C7> + 2/i, a 5 = -9^(<Ti - i><r 8 ) 
a 9 = 2j>(a 10 - i>a 8 ) + 2{d - ^z> 3 ) - 2hu 2 



g± = a 8 + h, 92 = -oio + 2i>a 8 + ^Cz> 2 + 2/ii>. 



(4.29) 



19 It may be verified that (|4.29|) is consistent with ( |4.8| ). In other words the entropy 
current determined by comparison with partition function agrees exactly with the non 
dissipative part of the entropy current determined from the requirement of positivity 
of divergence. 20 

5. CPT Invar iance 



In this section we explore the constraints imposed on the partition function Q2.17 ) and 



( p.26|) by the requirement of 4 dimensional CPT invariance. In Table 3 we list the 
action of CPT on various fields appearing in the partition function. 

• Parity even case: Using this table one easily see that, demanding CPT invari- 
ance of the action ( p. 17] ), the functions /1, f 2 , fz are even under CPT. Instead 
had we demanded only time reversal invariance, then the we would conclude that 

h = f2 = h = 0. 



19 The expression 2C\ was referred to as s 9 in 0. 

20 Note however that the entropy positivity method, in addition, determines two dissipative terms 
in the entropy current, and so, in that sense, carries more information about the entropy current. 
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Table 3: Action of CPT 



Parity odd case: Now demanding CPT invariance of the action Q2.26|) , we 
conclude that gi is odd function of A Q and hence it can not contain any constant. 
This in particular implies h = 0, since g\ = as + h. So the gauge non invariant 
piece in entropy current in ( |4.21|) vanishes once we demand CPT invariance. The 
function g 2 appearing in (|2.26|) is even function in A . It is also easy to see that 
the requirement of CPT invariance of the partition function forces C\ — 0. 



6. Discussion 

In this paper we have studied the equality type constraints between transport coeffi- 
cients for relativistic superfluids at first order in the derivative expansion. Our central 
result is that the constraints obtained from a local form of the second law of thermo- 
dynamics agree exactly with those obtained from a study of the equilibrium partition 
function. 

As the constraints obtained from both methods are numerous and rather involved 



in structure, the perfect agreement found in this paper strengthens the conjecture [[12 
that the constraints obtained from the partition function agree with those obtained 
from the local version of the second law of thermodynamics under all circumstances. 
It would be interesting to find either a proof for or a counterexample against this 
conjecture. 



In the special case that the superfluid is nondissipative, |19| has presented a frame- 
work for describing superfluid dynamics from an action formalism. It would be inter- 
esting to understand the connection of the formalism of [1T5] to that described in this 
paper. 

As we have explained above, a central object in our analysis was a local Euclidean 
action for the superconducting Goldstone field. In the neighborhood of a second order 
phase transition familiar Landau- Ginzburg action for the order parameter is the natural 
analogue of the Goldstone boson action utilized in this paper. It seems likely that 
the methods of the current paper generalize to the study of hydrodynamics in the 



neighborhood of second order phase transitions (see |^TJ for a review). It would be 
interesting to perform this generalization. 

Finally, in this paper we have discussed only the equality type constraints on 
nondissipative transport coefficients that follow from the local second law. We have 
neither discussed Onsager type equality constraints on dissipative coefficients nor the 
inequalities on dissipative coefficients. It is possible that these constraints follow the 
imposition of reasonable conditions (like stability) to time fluctuations about equilib- 
rium. We leave the study of time dependence to future work. 
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